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ABSTRACT 


Differential  energy  loss  (Mev/mg/cm1)  and  range  (mg/cm*)  data  have 
been  calculated  for  low  and  medium  energy  protons  and  alpha  particles 
in  the  following  eight  explosives:  RDX  (cyclotrimethylene  trinitramine), 
HMX  (cyclotetramethylene  tetranitramine),  TNT  (2,4,6-trinitrotoluene), 
PETN  (ptntaerythritol  tetranitrate),  tetryl  (2,4,6-trinitrophenyl  methyl 
nitramine),  lead  styphnate,  mercury  fulminate,  and  lead  azide.  The  well- 
known  Bethe  theory  was  followed  in  establishing  the  proton  ranges. 
Experimental  proton  differential  energy  loss  data  and  theoretical  compu¬ 
tations  were  used  to  establish  atomic  stopping  number  versus  energy 
curves  for  the  elements  hydrogen,  carbon,  nitrogen,  oxygen,  and  lead 
between  20  Kev  and  25  Mev.  From  these  curves,  molecular  stopping 
numbers  were  calculated.  The  alpha-particle  ranges  were  obtained  from 
the  established  proton  ranges. 
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INTRODUCTION 


In  the  course  of  investigations  of  heavy  particle  radiation  damage  :o 
explosives,  we  have  noted  the  absence  of  any  accurate  range-energy  rela¬ 
tions  for  these  materials  in  the  literature.  The  ranges  in  explosives  for 
alpha  particles  and  fission  fragments  were  calculated  in  a  Picatinny 
Arsenal  technical  report  (Ref  1);  however,  the  method  used  was  not  appli¬ 
cable  to  alpha  particles. 

Range-energy  relations  of  protons  in  explosives  can  be  calculated  using 
the  theory  of  Bethe  (Ref  2).  In  this  report,  Bethe’ s  method  was  used  to 
calculate  the  ranges  for  protons  of  energies  from  0.020  Mev  (for  the  sec¬ 
ondary  explosives)  or  0.050  Mev  (for  the  primary  explosives)  to  25  Mev  in 
eight  different  explosives  (See  Table  1,  p  6).  Alpha-particle  range-energy 
curves  were  then  derived  from  these  proton  ranges. 

DISCUSSION 

The  differential  energy  loss  (Mev/mg/cmJ)  for  protons  is  given  by: 

-i£=2»re4/-^\  !il(B/E)  (1) 

dR  \  me /  M 

Here,  e  is  the  electronic  charge  (for  these  calculations,  it  is  convenient  to  take 
e1  =  14.397  x  10- 14  Mev-cm);  Mp/me  is  the  ratio  of  proton  to  electron 
masses  (1836.6);  N0  is  Avogadro’s  number;  M  is  the  atomic  weight  of  the 
element  (or  molecular  weight  of  the  compound)  in  mg/mole;  B  is  a  dimen¬ 
sionless  quantity  called  the  atomic  (or  molecular)  stopping  number;  and  E 
is  the  energy  (Mev)  of  the  proton. 

By  rearranging  Equation  1,  the  atomic  stppping  number  for  each  element 
can  be  expressed  by 


B  =  -dE/dR  JL  E  (2) 

Mp  N„  2  ire4 

The  molecular  stopping  number  of  the  explosives  is  obtained  by  summing 
the  appropriate  atomic  stopping  numbers  (Refs  3  and  4)  as  follows: 

Bmol  <Xr  Y,  . >  =  'Batom  W  +  sBatom  00  + .  (3) 
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The  proton  ranges  Rp(E)  (mg/cm*)  were  calculated  by  graphical  integration 
of 


r"(E)-/e(tw  <4) 

Stopping  numbers  were  calculated  down  to  0.020  or  0.050  Mev  (roughly  the 
limit  of  the  available  data)  to  minimize  the  initial  energy  for  which  a  range 
would  have  to  be  estimated.  It  should  be  noted  that  even  rather  large  er¬ 
rors  in  the  stopping  number  over  this  energy  interval  do  not  appreciably 
affect  the  value  of  the  range  for  intermediate  energy  protons. 

The  ranges  at  0.020  Mev  or  0.050  Mev  for  these  explosives  were  esti¬ 
mated  by  assuming  (Ref  3)  that  the  form  of  -dE/dR  over  this  interval  is 
proportional  to  E  4. 

Hence,  integrating  Equation  4  yields: 

Rp  (E)  =  2Ke'/5  (5) 

The  values  of  K  were  determined  using  the  calculated  -dE/dR  values  at 
0.020  or  0.050  Mev  for  the  appropriate  explosive. 

Alpha-particle  ranges  were  then  obtained  from  the  established  proton 
ranges  using  (Ref  5) 


1.007  Ra  (3.972E)  =  Rp  (E)  +  8  mol  (6) 

Each  constant  (8moi)  represents  the  difference  in  range  caused  by  varia¬ 
tions  in  electron  capture  and  loss.  The  value  of  ^mo]for  each  explosive 

was  estimated  from  the  known  value  for  air  (Ref  5)  by  assuming  that  this 
effect  depended  only  upon  the  molecular  electron  density,  i.e.,  that  the 
product  of  5moi  and  the  electron  density  is  a  constant  for  all  materials. 

RESULTS 

Experimental  proton  differential  energy  loss  data  (-dE/dR)  are  available 
for  hydrogen,  carbon,  nitrogen,  oxygen,  and  lead  primarily  for  energies  of 
1  Mev  and  below.  These  data,  which  are  shown  in  Table  2  (p  7),  were 
used  to  establish  the  atomic  stopping  number  curves  for  hydrogen  and  car¬ 
bon,  nitrogen  and  oxygen,  and  lead  as  shown  in  Figures  1  through  3  (pp  12 
through  14),  respectively.  The  lower  limits  of  the  curves,  0.020  Mev  for  the 
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light  elements  and  0.050  Mev  for  the  lead,  were  chosen  arbitrarily  and  some 
minor  extrapolations  were  necessary  to  cover  this  range.  Whenever  the  ex¬ 
perimental  data  on  the  same  element  differed  significantly,  Whaling’s 
evaluation  (Ref  6)  was  followed.  These  curves  were  then  extended  to 
25  Mev  by  means  of  the  following  theoretical  computations: 

1.  The  calculations  of  Aron  and  Hoffman  (Ref  7)  were  used  for  car¬ 
bon,  nitrogen,  and  oxygen. 

2.  Hirschfelder  and  Magee’s  calculations  (Ref  3)  for  hydrogen,  which 
above  1  Mev  agree  well  with  Whaling’s  evaluation  of  the  experimental  data, 
were  used  here.  The  data  were  extrapolated  from  15  to  25  Mev  using  the  re¬ 
sults  given  by  Aron  and  Hoffman  for  hydrogen  as  a  guide  over  the  interval. 

3.  Sternheimer’ s  (Ref  8)  calculations  were  used  for  lead.  In  those 
cases  (hydrogen  and  carbon)  in  which  two  or  more  theoretical  curves  were 
available  over  this  energy  region,  the  ones  which  agreed  best  with  the  data 
of  Brolley  and  Ribe  (Ref  9)  for  the  4.4+  Mev  protons  wete  used.  Finally, 
where  necessary  (carbon,  nitrogen,  and  oxygen),  the  theoretical  data  were 
adjusted  by  a  constant  multiplicative  factor  to  agree  with  the  4.4+  Mev  data. 

The  molecular  stopping  numbers  of  protons  in  the  eight  explosives  were 
obtained  using  the  proton  atomic  stopping  number  curves  (Figs  1-3,  ppl2- 
14)  and  Equation  3.  These  molecular  stopping  numbers  are  given  in  Table  3 
(p  8).  The  differential  energy  loss  (-dE/dR)  for  protons  of  energies  up  to 
25  Mev  in  the  explosives  was  then  calculated  by  means  of  Equation  1.  The 
values  obtained  are  shown  in  Table  4  (p  9). 

Ranges  for  protons  of  energies  up  to  25  Mev,  calculated  by  graphical 
integration  of  Equation  4,  are  presented  in  Table  5  (p  10).  Figures  4,  5,  and 
6  (pp  15,  16,  and  17)  present  proton  range-energy  curves  for:  (1)  RDX  and 
HMX,  which  are  also  representative  of  the  other  secondary  high  explosives; 
(2)  lead  styphnate;  and  (3)  lead  azide,  which  is  also  representative  of 
mercury  fulminate.  The  data  for  mercury  fulminate  are  less  reliable  than 
the  rest  since  BHg  was  assumed  equal  to  Bpb.  Experimental  measurements 

(Refs  4  and  10)  indicate  that  the  additivity  of  atomic  stopping  numbers  does 
not  hold  for  protons  of  energy  less  than  0.150  Mev.  Although  ranges  are 
given  for  energies  below  0.150  Mev  in  Table  5,  both  the  effects  of  chemical 
binding  on  the  additivity  and  the  greater  experimental  error  in  this  region 
decrease  their  validity. 
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The  alpha-particle  ranges  for  energies  up  to  approximately  100  Mev,  ob¬ 
tained  from  the  established  proton  ranges  as  discussed  previously,  are 
presented  in  Table  6  (p  11).  Figures,  4,  5,  and  6  (pp  15,  16,  and  17)  also 
present  the  appropriate  alpha-particle  range-energy  curves. 
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TABLE  2 


Selected  Experimental  Proton  Differential  Energy  Lett  Data 
Presented  a*  Stepping  Numbers 


Proton  Energy, 
Mev 

Hydrogen0 

Stopping  Ngmbera 
Carbon**'*  Nitrogen"* 

for 

Oxygon? 

Load* 

.020 

.427 

— 

1.09 

.937 

— 

.030 

.700 

- 

1.82 

1.63 

— 

.030 

.732 

1.82 

2.02 

— 

— 

.040 

1.03 

— 

2.58 

2.33 

— 

.040 

1.05 

2.52 

2.86 

2.55 

— 

.050 

1.36 

- 

3.39 

3.08 

— 

.050 

1.34 

3.09 

3.72 

3.43 

— 

.060 

1.63 

— 

4.19 

3.79 

— 

.060 

1.62 

4.12 

4.57 

4.24 

— 

.070 

1.85 

— 

— 

4.48 

— 

.070 

1.86 

4.90 

5.49 

5.02 

— 

.075 

- 

- 

- 

- 

12.35 

.080 

1.87 

— 

— 

5.09 

— 

.080 

2.08 

5.51 

6.19 

5.74 

— 

.090 

2.27 

6.15 

6.87 

6.50 

— 

.100 

2.44 

6.79 

7.49 

7.19 

17.39 

.125 

- 

— 

— 

— 

22.47 

.150 

2.95 

9.16 

10.11 

10.13 

27.40 

.175 

- 

— 

— 

— 

32.12 

.200 

3.26 

10.62 

11.88 

12.30 

36.45 

.250 

3.48 

11.79 

13.08 

13.87 

43.58 

.300 

3.65 

12.79 

14.06 

15.05 

48.91 

.350 

3.80 

13.61 

14.84 

16.13 

53.41 

.400 

3.93 

14.28 

15.63 

17.13 

57.52 

.400 

- 

— 

— 

— 

57.53 

.450 

4.02 

14.94 

16.23 

17.80 

60.77 

.450 

- 

- 

— 

— 

61.52 

.500 

4.12 

15.43 

16.91 

18.50 

63.40 

.500 

— 

- 

— 

— 

64.80 

.550 

4.18 

15.98 

17.51 

19.29 

67.14 

.550 

- 

- 

— 

— 

68.06 

.600 

4.26 

16.43 

18.10 

19.86 

70.23 

.600 

- 

— 

— 

— 

70.74 

.650 

- 

- 

— 

_ 

73.91 

.700 

- 

- 

18.65 

— 

75.50 

.750 

- 

- 

— 

— 

78.70 

.800 

— 

- 

— 

— 

80.60 

.850 

- 

- 

— 

— 

83.15 

.900 

— 

— 

— 

— 

86.54 

.928 

— 

— 

21.28 

_ 

_ 

.950 

- 

- 

— 

— 

88.96 

1.00 

— 

- 

21.81 

_ 

91.97 

4.40 

- 

- 

— 

36.42 

— 

4.42 

— 

— 

32.52 

_ 

4.43 

— 

29.26 

_ 

_ 

_ 

4.44 

6.22 

- 

- 

- 

- 

*Hydrogen  and  oxygen  data  were  obtained  from  References  4,  9,  and  11. 

Carbon  data  were  obtained  from  References  4  and  9 

CThe  carbon  values  fo*  the  proton  energies  from  0.030  to  0.100  Mev  were  calculated  from 
the  data  and  method  given  in  Reference  4. 

^Nitrogen  data  were  obtained  from  References  4,  9,  11,  and  14. 

*Lead  data  were  obtained  from  References  12  and  13. 
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Fig  1  Stopping  numbers  of  protons  in  hydrogen  and  carbon 


13 


Fig  2  Stopping  numbers  of  protons  in  nitrogen  and  oxygen 
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Fig  6  Proton  and  alpha-particle  range-energy  curves  for  lead  azide 
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